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We consider the evolution of cavities within spherically symmetric relativistic fluids, under the 
assumption that proper radial distance between neighboring fluid elements remains constant during 
their evolution (purely areal evolution condition). The general formalism is deployed and solutions 
are presented. Some of them satisfy Darmois conditions whereas others present shells and must 
satisfy Israel conditions, on either one or both boundary surfaces. Prospective applications of these 
results to some astrophysical scenarios is suggested. 
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I. INTRODUCTION 

Many years ago, Skripkin |l| addressed the very inter- 
esting problem of the evolution of a spherically symmetric 
fluid distribution following a central explosion. As a re- 
sult of the conditions imposed by Skripkin a Minkowskian 
cavity should surround the centre of the fluid distribu- 
tion. 

Recently [_| , this problem was studied in detail, prov- 
ing that under Skripkin conditions (isotropic fluid with 
constant energy density distribution) the scalar expan- 
sion vanishes. It was further shown that the assumption 
of vanishing expansion scalar requires the existence of a 
cavity within the fluid distribution (of any kind). Next, it 
was shown in [3J] , that the Skripkin model is incompatible 
with Darmois junction conditions [_|. Also, the inhomo- 
geneous expansionfree dust models presented in Q are 
deprived of physical interest since they imply negative 
energy density distributions. 

For the reasons above, we turn in this paper to another 
kincmatical condition particularly suitable for describing 
the evolution of a fluid distribution with a cavity sur- 
rounding the centre. This consists in assuming the van- 
ishing of variation of proper radial distance between any 
two infinitesimally close fluid elements per unit of proper 
time. We shall explore here the consequences derived 
from this condition (hereafter referred to as the purely 
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areal evolution condition) . In particular we are interested 
in analytical models which even if are relatively simple 
to analyze, still contain some of the essential features of 
a realistic situation. It should be emphasized that we are 
not interested in the dynamics and the conditions of the 
creation of the cavity itself, but only in its evolution once 
it is already formed. 

We have two hypersurfaces delimiting the fluid. The 
external one separating the fluid distribution from 
a Schwarzschild or Vaidya spacetime (depending on 
whether we assume the evolution to be adiabatic or dissi- 
pative) and the internal one, delimiting the cavity within 
which we have Minkowski spacetime. It should be men- 
tioned that for cavities with sizes of the order of 20 
Mpc or smaller, the assumption of a spherically sym- 
metric spacetime outside the cavity is quite reasonable, 
since the observed universe cannot be considered homo- 
geneous on scales less than 150-300 Mpc. However for 
larger cavities it should be more appropriate to consider 
their embedding in an expanding Lemaitre-Friedmann- 
Robertson- Walker spacetime (for the specific case of void 
modeling in expanding universes see [_| , [_| and references 
therein). 

Thus, we have to consider junction conditions on both 
hypersurfaces. Depending on whether we impose Dar- 
mois conditions [|| or allow for the existence of thin shells 
0, different kind of models are obtained. In this paper 
we shall focus mainly on models satisfying Darmois con- 
ditions, although some models presenting thin shells will 
be briefly described too. 

For sake of generality we shall start our discussion by 
considering an anisotropic dissipative fluid (arguments 
to justify such kind of fluid distributions may found in 
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Q-pjl and references therein). A detailed description of 
this kind of distribution, as well as definitions of kinemat- 
ical and another important variables, are given in section 
II. The Darmois junction conditions on both, the inner 
and the outer boundary surface, are briefly discussed in 
section III. 

In order to understand better the physical meaning of 
the purely areal evolution condition, we shall discuss two 
different definitions of radial velocity of a fluid element, 
in terms of which both the expansion and the shear can 
be expressed, and which renders intelligible the origin of 
the term we use to denote such a condition in section IV. 

We shall next deploy all the equations required for cav- 
ity modeling under the purely areal evolution condition 
in section V. The specific case of cavities satisfying Dar- 
mois conditions, on both hypersurfaces, is treated in sec- 
tion VI, whereas cavities presenting shells are considered 
in section VII. 

Finally a brief summary of the results is presented in 
the last section and prospective applications of these re- 
sults are briefly mentioned. 



II. FLUID DISTRIBUTIONS AND 
KINEMATIC AL VARIABLES 

We consider a spherically symmetric distribution of 
fluid, bounded by a spherical surface Y,( e \ The fluid is 
assumed to be locally anisotropic, with principal stresses 
unequal, and undergoing dissipation in the form of heat 
flow (diffusion approximation). 

Choosing comoving coordinates inside H^ e \ the general 
interior metric can be written 

ds 2 _ = -A 2 dt 2 + B 2 dr 2 + R 2 (d9 2 + sin 2 6d<j) 2 ), (1) 

where A, B and R are functions of t and r and are as- 
sumed positive. We number the coordinates x° = t, 



and its shear <j a p by 



r, x 



9 and x 



Observe that A and B 



are dimensionless, whereas R has the same dimension as 
r. 

From (TTJ) we have that inside T,^ any spherical surface 
has its proper radius given by J Bdr and its areal radius 
byR. 

The matter energy- momentum T Q/3 inside Y,( e > has the 
form 

T a/3 = 0*+ P x)V a Vp+P±g a l3 + (P r -P±)x a Xl3+<laVp+V a qp, 

where fi is the energy density, P r the radial pressure, Pj_ 
the tangential pressure, q a the heat flux, V a the four- 
velocity of the fluid and \ a a unn - four- vector along the 
radial direction. These quantities satisfy 

V a V a = -l, V a q a =0, x a X a = h X a V a =0. (3) 

The four-acceleration a a and the expansion O of the fluid 
are given by 



(5) 



where 

h a /3 =g a p + V a V/3. (6) 

Since we assumed the metric ([T]) comoving then 

V a = A- x 5%, q a = qB^S?, X a = B' 1 ^, (7) 

where q is a function of t and r. From ^ with {7J we 
have the non zero component for the four-acceleration 
and its scalar, 



=-. a = (a«a a )i/ 2 = 4 E , (8) 



A[_ 

T 



AB 



and for the expansion 



1 IB R\ 



(9) 



where the prime stands for r differentiation and the dot 
stands for differentiation with respect to t. With (JT]) we 
obtain for the shear ([5]) its non zero components 

2 I 
o-ii = gB 2 cr, CT22 = eras sin" 2 9 = - -R 2 cr, (10) 

and its scalar 



where 



1 B R 
a= A\B-R 



(11) 



(12) 



Also observe that the shear tensor may be written as: 



(13) 



Sometimes it could be convenient [Tl|, [l2[ to express 
the energy-momentum tensor ([2]) in the form 

T ap = VVaVp + Ph al3 + n Q(3 + q (V a X0 + X a V fj ) (14) 
with 

1 



Pr + 2P| 

3 



IF" 9 = f hth^ ] - lh al3 h lS \ T lS = II | y Q y' 3 - lh al3 



X X 



a a = V a , p V^, e = v%, 



(4) 



n = p r -p 



Next, the mass function m(t,r) introduced by Misner 
and Sharp (see also [IH) reads 



R 3 D 23 R 
m = -R^ = - 



5 -(§• > 



(15) 



To study the dynamical properties of th e sy stem, let us 
introduce, following Misner and Sharp [13J, the proper 
time derivative Dt given by 



Dn 



l__d 

Adt' 



(16) 



Using (fir?)) we can define the velocity U of the collaps- 
ing fluid (for another definition of velocity see section IV) 
as the variation of the areal radius with respect to proper 
time, i.e., 

U = DtR < (in the case of collapse). (17) 
Then (fT5"]) . by using (fTT|). can be rewritten as 

1/2 



R' ( , 2m 

E = — = { 1 + U 

B \ R 



(18) 



and (|17p we obtain from (| 15[) 



m' = 4tt [fx + q—] R'R 2 , 



which implies 



m = Att 



R 2 R'dr, 



(19) 



(20) 



where we assumed a regular centre to the distribution, 
so m(0) = 0. 

It will be useful to introduce the Weyl tensor. Thus, 
let E a p denote the "electric" part of the Weyl tensor (in 
the spherically symmetric case the "magnetic" part of 
the Weyl tensor vanishes, H a p = 0) defined by 



E a /3 = CauPvV^V" , 

which may be written as 



E a /3 = £ [ XaX/3 ~ \h a p ) , 



Using field equations (see [2j for details) ) with (fTo]) where 



(21) 



(22) 
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R" 




R'\ 
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2~R2' 



(23) 



Using field equations (see for details) and the defi- 
nition of mass function (|15l) we may write £ as 

£ = M^-Pr + P±)-^- (24) 



III. THE EXTERIOR SPACETIME AND 
JUNCTION CONDITIONS 



where M(v) denotes the total mass, and v is the retarded 
time. The matching of the non-adiabatic sphere to the 
Vaidya spacetime, on the surface r = r £ ( e ) = constant, in 
the absence of thin shells, where Darmois conditions hold, 
is discussed in 9, 15 17]. This requires the continuity of 
the first and the second fundamental forms through the 
matching hypersurface, producing 



Outside £( e ) we assume we have the Vaidya spacetime 
(or Schwarzschild in the dissipationless case), i.e., we as- 
sume all outgoing radiation is massless, described by 



1 



2M(v) 



dv — 2drdv 



l (d8 2 



sin 



2 ), 
(25) 



m(t,r) S i e> M(v), 



(26) 
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and 



or, by using (0 and ([T2j) , 



Q = 



(28) 



where = means that both sides of the equation are eval- 
uated on E^ 6 ) and L denotes the total luminosity of the 
sphere as measured on its surface and is given by 



L = Lnc I 1 — 



2m 



and where 



r dv 



dM 
dv 



(29) 



(30) 



is the total luminosity measured by an observer at rest 
at infinity. 

From (|27|) and field equations one obtains 



q = P r 



(31) 



In the case when a cavity forms, then we also have to 
match the solution to the Minkowsky spacetime on the 
boundary surface delimiting the cavity. If we call £W 
the boundary surface between the cavity and the fluid, 
then the matching of the Minkowski spacetime within the 
cavity to the fluid distribution, implies 



m(t,r) = 0, 



E<« p 
q = P r 



(32) 



(33) 



However, since we are assuming our cavity to be empty, 

S (i) 

then L = 0, which implies 



y(i) y(i) 

q = P r = 0. 



(34) 



If we allow for the presence of thin shells on £W and/or 
£( e ' , then we have to relax the above conditions and allow 
for discontinuities in the mass function Q. 



IV. TWO DEFINITIONS OF RADIAL 
VELOCITY AND THE PURELY AREAL 
EVOLUTION CONDITION 

In section II we introduced the variable U which, as 
mentioned before, measures the variation of the areal ra- 
dius R per unit proper time. 

Another possible definition of "velocity" may be intro- 
duced, as the variation of the infinitesimal proper radial 
distance between two neighboring points (51) per unit of 
proper time, i.e. Dt(51). Then, it can be shown that 
(see 0] for details) 



D T (5l) 1 , 



(35) 



D T (5l) 
51 



B 

AB' 



Then with (|§J}, (fl2]). ([P?]) and jSSJ we can write 



D T (6l) D T R _ D T (5l) 



51 



R 



51 



U 
R : 



(36) 



(37) 



and 



e 



D T (5l) 2D T R D T (5l) 



51 



R 



51 



Thus the "circumferential" (or "areal") velocity U, is re- 
lated to the change of areal radius R of a layer of matter, 
whereas Dt(51), has also the meaning of "velocity", be- 
ing the relative velocity between neighboring layers of 
matter, and can be in general different from U. 

In 21 it was shown that the condition O = requires 
the existence of a cavity surrounding the centre of the 
fluid distribution. 

Let us now consider the condition Dt(51) — 0, but 
U 7^ 0. From the comments above it is evident why we 
shall refer to it as the purely areal evolution condition. 

Now, if D T (5l) = then it follows from J37J and ([58]) 
that © = — 2a, feeding this back into the (01) component 
of the Einstein field equations (see [2j) for details) we get 



<tR' 
~R 



4irqR' 



whose integration with respect to r yields 

C(i) 4tt r rr' 



R 



R 



o 



dr. 



(39) 



(40) 



where £ is an integration function of t. It should be ob- 
served that in the case where the fluid fill all the sphere, 
including the centre (r — 0), we should impose the reg- 
ularity condition £ = 0. However since we consider the 
possibility of a cavity surrounding the centre, such a con- 
dition is not required. 

On the other hand, (|4H|) with (|3"T|) implies 



U=-( 



4:r / q — dr. 
lo 



E 



(41) 



Thus, in the non-dissipative case the purely areal evo- 
lution condition implies that U = U(t). This condition 
is clearly incompatible with a regular symmetry centre, 
unless U = 0. Therefore if we want the purely areal evo- 
lution condition to be compatible with a time dependent 
situation (U ^ 0) we must assume that either 

• the fluid has no symmetry centre, 
or 

• the centre is surrounded by a compact spherical 
section of another spacetime, suitably matched to 
the rest of the fluid. 
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Here we shall discard the first possibility since we 
are particularly interested in describing localized objects 
without the unusual topology of a spherical fluid without 
a center. Also, within the context of the second alterna- 
tive we have chosen an inner vacuum Minkowski spherical 
vacuole. 

Let us now consider the dissipative case. Then assum- 
ing the purely areal evolution condition, if the fluid fills 
the whole sphere (no cavity surrounding the centre) , and 
we have a symmetry centre, we have to put £ = and 
(|4Tj) becomes 



that the condition B — requires for consistency, the 
existence of a cavity surrounding the centre. 
Then the Einstein field equations become 



R 



1 R\ R" ( R'\ 1 



2 R' RA' 
87Tq =A\R~RA 



(46) 



( = -I t / q dr, 
Jo & 



(42) 



which is not incompatible with a regular symmetry cen- 
tre. In this case we shall assume in an ad hoc manner that 
a cavity surrounds the centre. However, this assumption 
is somehow suggested by the following qualitative argu- 
ment. 

In the case of an outwardly directed flux vector (q > 0) , 
all terms within the integral are positive and we obtain 
from §E§ and (02) that 9 > and U > 0. Now, dur- 
ing the Kelvin- Hclmholtz phase of evolution [18j . when 
all the dissipated energy comes from the gravitational 
energy, we should expect a contraction, not expansion, 
to be associated with an outgoing dissipative flux. In- 
versely, an inwardly directed flux (q < 0) (during that 
phase) would produce an overall expansion, not a con- 
traction as it follows from (|42"T) . 

Thus we have seen that the purely areal evolution con- 
dition appears to be particularly suitable to describe the 
evolution of a fluid distribution with a cavity surrounding 
the centre. 

Finally observe that using (|57|) and ([55)1 in (TT5|) . the 
purely areal evolution condition can be expressed in a 
covariant form as: 



9 ( h 



(43) 



8nP r 



1 

"A* 



'r 



} A R\ R 

j a~r)r 



A' R'\ R' I 
2 A + R )r-&> (47) 



1 R AR\ A" R" A' R' 
^ = -TAr-AR ) + ^ + li + AR^ 



and the non trivial components of the Bianchi identities, 



T. 



-a/3 



0, become, 



f 

A 



A + 2( M + P ± )- 



q' + 2q 



(AR)' 
AR 



0, (49) 



Rf 



-2(P r -P ± )-=0. 



(50) 



In the particular geodesic case we have A' = — > A = 
1. Then the field equations (|45ti48|) read 



In the next section we shall consider some models. 



V. MODELS OF CAVITIES 

We shall now study the general properties of models 
satisfying the purely areal evolution condition. 

The general picture is similar to that proposed by 
Skripkin [T[, namely, an explosion at the center initiates 
an overall expansion of the fluid, creating a cavity sur- 
rounding the centre. The difference here is that we shall 
not assume 9 = but instead, Dt(51) = 0. Thus we 
have because of (|36|) B = (but R ^ ) which means 
that B = B(r) and it can be chosen 



B = 1, 



(44) 



with no loss of generality. As mentioned before, the phys- 
ical appeal of this kind of models stems from the fact 



87T/1 



R 
R 



8nP r 



and it follows with (IT5 



- 2- 



R" (R! 



R 



R 



1 




2nn=jp+ 2n(P r - 2P±), 



or, using 



2tt^ = 27r(P r - 4P X ) - £■ 



(51) 



8nq = 2^, (52) 



w (53) 



8nP ± = -- + —, (54) 



(55) 



(56) 
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Hence, from ([55)1 and (jSij)) we have for conformally flat 
spacetime £ — 0, and geodesic fluids with isotropic pres- 
sures P r = P± = P, 



R 3 



4ttP = 0, 



(57) 



implying that such models, if they satisfy Darmois con- 
ditions, must be dissipative (otherwise M = 0), and ab- 
sorbing energy (q < 0, otherwise m < 0). 



VI. MODELS SATISFYING DARMOIS 
CONDITIONS 

Let us now consider some simple cases in order to find 
analytical models which do not present thin shells on 
either E^ e ' or £w, but holding Darmois conditions. 

The simplest models of this kind we have found are 
non-dissipative. Thus, considering q — then from (|46[) 
after integration we have 



A — — 



(58) 



where hi(t) is an arbitrary function of t. Reparametriz- 
ing t, we may choose without loss of generality 



hi — i? S (i) , 
which amounts to choose 

Asm = 1. 



(59) 



(60) 



Observe that for all these models the velocity U is the 
same for all fluid elements, since, as it follows from (fT7| 
and JS 



U = h 1 =R^ ) . 



(61) 



This fact was already brought out in the previous section 
from (j4Tj) . 

Substituting (J58J) into (|4"5]). (H7|) and 0S} we obtain, 
using ([16]), dSU and ([59]) 



87rF r 



1 

P 2 ' 
1 



8 7 r^ = - 1 ^(2i? J R" + i?' 2 -i?| (l) -1), (62) 



■Dt^' 2 -^,,-!)], (63) 



1 



8?rPj - = T^^ d t£RR" + R' 2 - Rim - !)■ (64) 

From dm and d® it is clear that 

DrifiR 2 ) 



Pi = -- 



2RR 



(65) 



Calculating the mass function (IT51) with dS]) and 
it yields 



m = -^(i?' 2 -7?| w -l), 



(66) 



which allows to reexpress (|6"3"|) like 

TT7. 

47rP r 



i? 2 i? 



(67) 



The boundary conditions P r = and P r 



automatically satisfied with m 

m = 0. 

With dMl we obtain for (EI 



are 

M = constant and 



£ = 



R 



■Dn 



1 

P 2 



(2RR" - R' 2 + 7? 2 (!) + 1 



(68) 



4i? sw 

We shall next specialize to some particular cases. 

A. Conformally flat models 

If we assume the spacetime between r = r S ( e > and 
r = r S (i> to be conformally flat £ = 0, then it follows 
from UHl) that 

2PP" + R' 2 - AP 2 + i?| (<) +1 = 0, (69) 

where /i(r) is an arbitrary function of r, and integrating 
again we obtain 



R u = R 



hdR + h- t 



R^m + !) 



(70) 



where /i2(t) is an arbitrary function of t. Comparing ([66 
with ([70|) it follows that 



m = -^-( / ./;,</./? //, j . 



(71) 



therefore /12 may be obtained from the junction condition 
(|32|) . producing 



h 2 = - fidR. 



(72) 



Thus all models of this kind are defined by a single func- 
tion fi(r) which should be chosen such as to satisfy the 
remaining Darmois conditions. 

Before proceeding further, the following remark is in 
order: all spherically symmetric conformally flat space- 
times (without dissipation) and isotropic fluids are shear- 
free (see eq.(78) in [l(|), but this is no longer true for 
anisotropic fluids (see eq.(84) in [I3|). Therefore the 
models to be considered here are necessarily anisotropic. 



B. Models with vanishing tangential stresses 



Assuming P±_ = we have from (|65[) after integration 



R 2 ' 



(73) 
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where /2(f) is an arbitrary function of r. Substituting 
(|73| into dSU it yields 



2RR" + R' 2 + 8tt/ 2 - R 
and with (po]) it becomes 

to' = 47r/2-R' 



2 

£(<) 



1 = 0. 



(74) 



(75) 



In order to obtain models we have to assume specific 
form of the mass function or the energy density distribu- 
tion. As an example let us assume 



/2 = C\ = constant > 0, 
producing because (|20|) and (|73| 
to = 4irci(R - 

and implying 

M = 47rci(i? s( e) -R sW ), 

i? S (s) = Rj:(i), A S (<) = A S (e) = 1. 



(76) 



(77) 



(78) 
(79) 



Also, from (|17]), (J73]) and (J771) it follows 



= M Z !)• 

-ft 



Next, substituting (1771) into (1551) we have 



where 



a(t) = + 1 - 87TCi, /3(*) = 87rcii? s(l ) 



(80) 



(81) 



(82) 



and after integration 

[ai?(ai?+/3)] 1/2 -^ln[(ai?) 1 / 2 + (ai?+/?) 1 / 2 ] = a 3 / 2 [r-r ], 

(83) 

where ro(t) is an arbitrary function of t. Evaluating (|53")) 
on £w we obtain 



a(i?| (!) + 1) 



1/2 



87TCi In 



« 1/2 + (4<o + 1) 1/2 



4ttci In i? S> R a 3/2 (r-r ). 



(84) 



This is a first order differential equation for which 
may be solved for any function ro(t). The result of this 
integration, together with (1831) provide all the informa- 
tion required to obtain the t and r dependence of all 
physical and metric variables. Observe that the energy 
density is always positive and regular everywhere within 
the fluid distribution. Also, choosing ro(t) such that 

< ^ — 1 < 1 we assure that the pressure is posi- 
tive and smaller that the energy density. 



VII. MODELS WITH THIN SHELLS 



As it is apparent from the discussion above, the ful- 
filment of Darmois conditions on, both, £( e ) and £W 
severely restrict the possible models of cavities. There- 
fore it might be pertinent to relax Darmois conditions 
and work within the thin wall approximation, which al- 
lows for the existence of discontinuities of the mass func- 
tion across X^ 6 - 1 and/or X^ (for models of voids within 
the thin wall approximation see [l9|-(2l1| and references 
therein). 



A. Non-dissipative geodesic model 

The simplest model of this kind (under the purely areal 
evolution condition) corresponds to a Lemaitre-Tolman- 
Bondi (LTB) spacetime (22|- [HI, whose general line ele- 
ment is given by 

D/2 

ds 2 _ = -A 2 dt 2 + dr 2 + R 2 {d6 2 + sin 2 6d(j) 2 ), (85) 

1 — K 

where K is a function of r. Then, the purely areal evolu- 
tion condition applied to (|85|) produces 

i?' = (l-iq 1/2 , (86) 

and since all LTB spacetimes are geodesic and dissipa- 
tionless, (|52|) implies 

R = h 3 +f 3 , (87) 

where h^{t) and fs(r) are arbitrary functions of t and r 
respectively. Then for the mass function (p~5]l we obtain 

m = | (hl-ft + l). (88) 

Imposing Darmois conditions on £W we obtain from 
([32]) with (18611881) . 

/i 3 = (-X sW ) 1 / 2 < + c 2 , (89) 
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where C2 is an arbitrary constant and we have to assume 
the condition < 0. Once we have imposed Darmois 

conditions on E^- 1 , it follows that we have to assume the 
presence of a shell on £( e ). Indeed, the continuity of the 
mass function m on £( e ) implies from (1261) 



R 



2M 



hi-fi 2 



(90) 



which, because of and the fact that M is constant, 
produces 



R = 0, 



implying, because of ([87[) . 



0^i? = 0, 



(91) 



(92) 



thereby invalidating the continuity of mass function 
across Y,( e \ Thus all models of this kind should ad- 
mit a thin shell on £( e ). Alternatively, we may assume 
discontinuities of the mass function on £w or on both 
boundaries. 



B. Geodesic radiating dust models 

In an increasing order of complexity, let us consider the 
next simplest possible situation, namely, geodesic a — 0, 
dust P r = P± = with dissipation q ^ 0. It should 
be stressed that in the dissipative case, the pure dust 
condition, P r = P±_ = 0, does not imply vanishing four- 
acceleration a, as it can be seen from (E 

From (l54"l) it follows at once 



R — i?", 

whose general solution is of the form 

R = c 3 ty(t + r) + c A <S>{t - r), 



(93) 



(94) 



where C3 and C4 are constants. Considering (|93l) together 
with (j5"Tj) and (f53"| produces 



2-7T/1 



R" 



R 



R R 1 
and (j50|) in our case reduces to, 

q + 2q- = 0, 

implying 

u 



(95) 



(96) 



(97) 



where fi(r) is an arbitrary function of r. 

We shall now assume that our model satisfies Darmois 
conditions on then from (l34l) it follows 



s (0 

U = 0. 



(98) 



Next, from ©, ([171 and ppl) we have 



(99) 



R ~-Rj' 



and evaluating it at the cavity boundary by using ([32 
we obtain 



R = c 5 t + c 6 



(100) 



where C5 and cq are constants. Observe that this is con- 
sistent with ([Ml). Also it follows from J95l) that 



At = 0. 



(101) 



A further restriction on f± may be obtained from (|49p . 
which in our case reads 



R R 1 
/i + 2/i- +q' + 2q— = 0. 



(102) 



Evaluating (I102p on the boundary of the cavity and using 
(j98l) and (fT0T|) . we have 



/i E " 0. 



(103) 



Once we have assumed that Darmois conditions are 
satisfied on S™, then it follows that they are violated 
on £( e ). Indeed, if we assume continuity of the mass 
function across E^ e \ then evaluating (|99[) on £( e ) and 
using (|2"6"|) we obtain 



- £<<=> M 

R - ~W 
which can be integrated to obtain 



• S Ce) (2M 

R = ± U +c7 



1/2 



(104) 



(105) 



where C7 is a constant. Combining (|97[) and (|102[) we 
obtain 



(106) 



/ii? 2 = / B - tfi 



where /s(r) is an arbitrary function of r. Then using (|95[) 
and (fTUtl) in (fT06]) and evaluating at £( e ' it follows that 



71/ 



Mfs-tfi)' 

and by using p05p becomes 



i? = constant, 



(107) 



(108) 



thereby ruling out the possibility of the continuity of the 
mass function across T,( e \ Observe that in this particular 
case, there should be always a shell on T,^ . 

The only constraints imposed on functions fi, \1/ and 
$ are, that they must be regular and positive, so that the 
regularity and positiveness of /i is assured for all values 
of t and r. 
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C. Non geodesic models 

So far all presented models within the thin wall ap- 
proximation have been geodesic, therefore it would be 
instructive to present a non geodesic model. For that 
purpose, we shall invoke an ansatz which proved to be 
very useful for describing dissipative collapse (l6| . 

Thus, let us assume 



A' 



;2 
'A 



R' Q \ R' 



Kq J Ho K 



1 



8nP ±0 = 



A" Ft" A' R' 



A R A Ro 
With (|109H112I) we can rewrite (|45H48p like 



(111) 



(112) 



A = A (r), R = R Q (r)g(t), 



(109) 



where we take Aq and Ro to describe a static anisotropic 
perfect fluid whose energy density no and anisotropic 
pressures P r o and P±o are given by 



87T/! 



, R o 

'Ro 



Ro 



1 

Hq 



(110) 



87T/i = KfJLQ 



8nq 



1 

Al 



o V5 



-(- 



2 A' \ g 



A \R A a J g 



(113) 



(114) 



87rP r 



2 9 - 



— ( —- 



(115) 



8irP± = 8wP ±0 



i g 



A> 



6 9 



(116) 



All models of this kind present shells in either E^-* 
or £( e ). Indeed, evaluating (I115|) on S™ and assuming 
Darmois conditions there, we obtain 



2gg + g 2 - c 2 (g 2 - 1) = 0, 



where 



co 



Rq 



Then integration of T[) produces 



• 2 2 g , 

5 = Cn I — - 1 
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(117) 



(118) 



(119) 



where cs is a constant. Finally, evaluating the mass func- 
tion (JTSJ on EW, and using (I119p it follows that 



g(t) = constant. 



(120) 



implying the necessary violation of Darmois conditions 
on E« . 

Choosing a physical meaningful static ("seed") solu- 
tion, it is not difficult to choose function g such that 
standard energy conditions are satisfied. Thus for exam- 
ple if we demand g 2 < 1 and - > 0, we assure those 
conditions for the models. 



VIII. CONCLUSIONS 

We have studied in detail the consequences emerging 
from the purely areal evolution condition. It has been 



shown that such a condition is particularly suitable for 
describing the evolution of a fluid distribution endowed 
with a cavity surrounding the centre. 

All equations governing the dynamics under such con- 
dition have been written down and some models have 
been presented. 

Some of them satisfy Darmois conditions on both hy- 
persurfaces, £W and E^ e )), precluding thereby the ap- 
pearance of shells on either of these hypersurfaces. More 
simple models result from relaxing Darmois conditions 
and adopting Israel junction conditions across shells. 

One possible application of the presented results which 
comes to our minds, is the modeling of evolution of cos- 
mic voids. Indeed, the cavity associated to the purely 
areal evolution condition might be considered as a void 
precursor. 

Voids are, roughly speaking, underdensity regions in 
the larg e-scale matter distribution in the universe (see 
[25|-|33|| and references therein). 

Their relevance stems from the fact that, as stressed 
in [34j, the actual universe has a spongclike structure, 
dominated by voids. Indeed, observations suggest [35[ 
that some 40-50% of the present volume of the universe 
is in voids of a characteristic scale 30 h _1 Mpc, where 
h is the dimensionless Hubble parameter, Hq=100 h km 
sec -1 Mpc . However voids of very different scales may 
be found, from minivoids (36j to supervoids (37j | . 

It should be emphasized that in general voids are nei- 
ther empty nor spherical, either in simulations or in deep 
redshift surveys. However, for simplicity they are usually 
described as vacuum spherical cavities surrounded by a 
fluid. 

However we are aware of the fact that cold dark matter 
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at scales of the order of tens of Mpc is non-collisional, so 
that pressure and heat flux terms are negligible. There- 
fore, excluding the LTB case, it is not likely that our 
solutions could be used as toy models for cosmic voids. 

Possibly our solutions could be applied as toy models of 
localized systems such as supernova explosion models. It 
is worth mentioning that for these scenarios, the Kelvin- 
Helmholtz phase is of the greatest relevance (38l |. 

At any rate our purpose here has not been to generate 
specific models of any observed void, but rather to call 
the attention to the potential of the purely areal evolution 



condition for such a modeling, providing all necessary 
equations for their description. 



Acknowledgments. 

We would like to thank an anonymous referee for his 
very thorough report, which substantially helped us to 
improve our manuscript. This research has made use of 
NASA's Astrophysics Data System. 



V. A. Skripkin Soviet Physics-Doklady 135, 1183 (1960). [19 

L. Herrera, N. O. Santos and A. Wang Phys. Rev. D 78, [20 

084026 (2008). [21 

L. Herrera, G. Le Denmat and N. O. Santos Phys. Rev. [22 

D 79, 087505 (2009). [23 

G. Darmois Memorial des Sciences Mathematiques [24 

(Gauthier-Villars, Paris, 1927) Fasc. 25. [25 

W. Bonnor and A. Chamorro Astrophys. J. 361, 21 [26 

(1990). [27 

R. Torres Class. Quantum Grav. 22, 4335 (2005). [28 
W. Israel II Nuovo Cimento B 56, 1 (1966). 

L. Herrera and N. O. Santos Phys. Rep. 286, 53 (1997). [29 
L. Herrera and N. O. Santos Phys. Rev. D 70, 084004 

(2004). [30 

L. Herrera, A. Di Prisco, J. Martin, J. Ospino, N. O. [31 

Santos and O. Troconis, Phys. Rev. D 69, 084026 (2004). [32 

R. Maartens, astro-ph/ 9609119. [33 

L. Herrera, A. Di Prisco and J. Ospino [34 

arXw:1001.3020vl. [35 
C. Misner and D. Sharp Phys. Rev. 136, B571 (1964). 

Cahill M. and McVittie G. J. Math. Phys. 11, 1382 [36 
(1970). 

N. O. Santos Mon. Not. R. Astron. Soc. 216, 403 (1985). [37 
W. B. Bonnor, A. Oliveira and N. O. Santos, Phys. Rep. 

181, 269 (1989). [38 
R. Chan Mon. Not. R. Astron. Soc. 316, 588 (2000). 
R Kippenhahn, A. Weigert Stellar Structure and Evolu- 
tion, Springer Verlag, Berlin, (1990). 



K. Lake and R. Pirn Astrophys. J. 298, 439 (1985). 
R. Pim and K. Lake Astrophys. J. 304, 75 (1986). 
R. Pim and K. Lake Astrophys. J. 330, 625 (1988). 

G. Lemaitre Ann. Soc. Sci. Bruxelles A 53, 51 (1933). 
R. C. Tolman Proc. Natl. Acad Sci 20, 169 (1934). 

H. Bondi Mon. Not. R. Astron. Soc. 107, 410 (1947). 
Ya. Zeldovich Astron. and Astrophys. 5, 84 (1970). 

P. J. E. Peebles arXiv.astro-phj '0101127. 

P. J. E. Peebles Astrophys. J. 557, 495 (2001). 

F. Hoyle and M.S. Vogeley Astrophys. J. 566, 641 

(2002). 

M. Serpico, R. d'Abrusto, G. Longo and C. Stomaiolo 

Gen. Rel. Grav. 39, 1551 (2007). 

R. Giovanelli arXw:0909.3U8vl. 

A. Odrzywolek Phys. Rev.D 80,103515 (2009). 

J. Gaite arXiv:0911.1871vl. 

R. van Weygaert and E. Platen arXiv. 0912.2997. 
D. L. Wiltshire arXiv:0712. 39841. 

F. Hoyle and M. S. Vogeley Astrophys. J. 607, 751 
(2004). 

A. V. Tikhonov and I. D. Karachentsev Astrophys. J. 
653, 969 (2006). 

L. Rudnick, S. Brown and L. R. Williams 
arXiv.0704.0908. 

A. Burrows and J. Lattimer Astrophys. J. 307, 178 
(1986). 



